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METHOD OF LEAST SQUARES APPLIED TO ESTI- 
MATING ERRORS IN COAL ANALYSIS. 


By J. D. Davis and J. G. Farrcnip. 


INTRODUCTION. 


In this paper, which is published by the Bureau of Mines as a 
contribution to the increase of efficiency in the utilization of fuel, 
the authors have sought to describe a definite method for estimating 
the probable limits of error in analyzing samples of coal, and to show 
how this method, that of least squares, should be applied. 

Bailey * has shown that the method of least squares is applicable 
to the problem of sampling coal; it can be used to determine the 
limits of accuracy of sampling and for preparing standardized 
methods of sampling, but so far as the authors know the accuracy 
of the ordinary methods of analyzing coal has not yet been investi- 
gated by a mathematical method. It is customary to estimate the 
limits of accuracy of analysis by observing the agreement of duplicate 
determinations; that is, from general experience with a given kind 
of analysis. Hence, as some sets of duplicates are found to agree 
more closely than others, this question arises: What is the widest 
variation and what is about the average variation reasonably to be 
expected? Or, in other words the question is: What is the magnitude 
of the “probable error” of a single determination and what is the 
Maximum error to be expected, say, once in 1,000 determinations ? 

Although by the method of least squares direct expressions may 
be obtained for the magnitude of the errors in the work studied, the 
authors do not wish to imply that ideas formed from long experience 
are of no value; the theory of probability and the method of least 
squares are based on just such experience. Leplace has said “the 
theory of probability is nothing but common sense reduced to calcu- 
lation.”” The chief value of the theory is that it furnishes a ready 
means by which the individual may obtain for himself a concrete 
expression of the precision of his own work and at the same time a 
convenient form of statement for the benefit of others. 

The application of the theory of probability to obtain expressions 
for the accuracy of the proximate analysis of coal is complicated 


@ Bailey, E.G. Accuracy in sampling coal: Jour. Ind. and Eng. Chem., vol. 1, March, 1909, p. 161. 
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through the fact that the same methods are used for all kinds of coal 
and for no two coals will the several constituents, as stated, represent 
exactly the same quantities chemically. Consequently the error is 
often afunction of the stated chemical composition of the constituent. 
For example, if the mineral matter of a coal contains carbonates in 
varying percentages the probable error of the ash determination 
will be greater if the percentages of the carbonates are not determined 
and the proper correction made.* 

Further, it is obvious that the analytical error will be to some 
degree a function of the quantity of the constituent determined. A 
rigid treatment would involve the examination of all classes of coals, 
and of coals from every bed represented, and would require an enor- 
mous amount of analytical work. 

The authors have attacked the problem by first carefully working 
out the limits for a few coals that are selected as typical and then 
comparing the results so obtained with data compiled from numerous 
check determinations taken consecutively from the laboratory 
records. The number of check determinations was usually 1,000. 
The experimental work was all done by the same analyst; hence the 
effect of the personal equation was eliminated from the results. On 
the other hand the data taken from the records include the work of 
several analysts, so that the personal equation was here involved. 
Further, the experimental work for a given result took only a day or 
two, whereas the data represented in the compilation from the 
records represented at least a month’s analyses, so that variations 
chargeable to changing conditions in the laboratory probably were 
involved in the recorded results used. 

The committee on coal analysis of the American Society for Testing 
Materials has designated ° limits of accuracy for each determination 
involved in coal analysis, and as this information represents the 
consensus of opinion of leading American chemists, the limits are 
given to serve as a standard of reference for the work that follows. 
The accuracy of the method of preparing laboratory samples should 
be checked frequently by resampling the rejected parts of the large 
(3 to 5 pounds) sample and preparing duplicate laboratory samples. 
Results for ash in duplicate samples should agree within the following 


limits: 
Limits of error of ash determinations for duplicate samples. 
Per cent. 
No Carhonaten: Presenteas ase Seiesien wry swag ka Smetnivelne esas sore lease 0.4 
Considerable carbonates and pyrite present............2.-..2--.--- ae 


More than 12 per cent ash aud containing considerable carbonate. 


Tune, 1913, p. 523, 
6 American Society for Testing Materials, Yearbook for 1915, p. 596. 
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The ‘allowable variations” in analysis of small laboratory samples 
are as follows: 


“ Allowable variations’’ in coal analysis, as designated by American Society for Testing 


Materials. 
Same Different 
Moisture determination. analyst, analysts, 
i per cent. per cent. 
Moisture under 5 percent..............22220002- eee cece ee 0.2 0.3 
Moisture over 5 per cent... .....-..2.. cece ceed cece eee eee 0.3 0.5 
Volatile matter determination. 
Bitiiminowecoale: <5 862 522 cG8<.0.4 sont Dawes dine ss Sloane's 0.5 1.0 
Meignitegssisseas s2i035 233 oh edsitonecaaustedatvescosteconte 1.0 2.0 
Ash determination. 
No carbonates present...........-...-.20-2+-+e005- See sea 0.2 0.3 
Carbonadtes presenit, 10.32 2225 ose camila ace cgusicadiwpesicaas 0.3 0.5 
More than 12 per cent ash and containing carbonates and 
PYM On s2 esos SVS cis neta eces Seearce teesoes wlan wetes 0.5 1.0 
Sulphur determination. 
Sulphur under 2 percent. ..223 <scdsccescacsesessea veins dare 0.05 0.10 
Sulphur over 2 per cent « .~2- -...siessebcicecs coeticcsscecues 0.10 0. 20 
Calorimeter determination. 
AML Cogley ste casctheted dedsceesatniee Pennies ener 0.3 0.5 


The term ‘‘allowable variations” should be construed not strictly 
in the sense that agreement within the prescribed limits should always 
be obtained with ordinary care in manipulation, but rather that the 
great majority of check determinations will agree within the limits 
given; the limit itself should be taken merely as a criterion of the 
desirability of repeating the determination. 


PREPARATION OF THE SAMPLE OF COAL FOR ANALYSIS. 


Samples of coal for analysis are received by the Bureau of Mines 
in one or more galvanized-iron cans containing when full about 3 
pounds each. The method of reducing a sample to a size and bulk 
convenient for making the analysis is as follows:4 

The entire sample, the coarsest particles of which will pass a 4-mesh 
screen (,;-inch holes), is first crushed by a roll crusher, a grab sample 
for the total moisture determination is taken, and the whole is put 
through a 20-mesh sieve. Any particles too coarse to pass the sieve 
are again run through the crusher. The bulk of this 20-mesh sample 
is now reduced to about 6 ounces by a Taylor & Brunton splitter 
(Pl. I). The 6-ounce sample is further crushed by a planetary disk 
crusher after which it is divided on the splitter. The resultant 
sample (about 3 ounces), which will all pass a 60-mesh sieve, is placed 
in a 4-ounce rubber-stoppered bottle for analysis. 


@Stanton, F.M., and Fieldner, A.C., Methods of analyzing coal and coke: Tech. Paper 8, Bureau of 
Mines, 1913, p. 1; Pope, G.s. , Sampling coal deliveries: Bull. 63, Bureau of Mines, 1913, p. 33, 
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As the precision attainable by this method of reduction depends 
largely on the distribution of the impurities in the sample, the authors 
sought at the outset to gain some idea of the variation of impurities 
in different parts of the sample. The separation of impurities by 
the gravity or washing method was considered, but as that method 
seemed particularly adapted to the somewhat complete separation 
of extraneous ash (slate, pyrite, or bone), the variations indicated 
would hardly be approached in the method of reduction applied. It 
seemed reasonable to assume that some of the materials making up 
the sample would crush more easily than others; that the coal sub- 
stance itself would perhaps crush more easily than the impurities. 
If the ash determination be used as a measure of the impurities 
present in different parts of a sample, the maximum variation 
attainablo by sizing would then be evident, and would serve as an 
indication of the difficulty of obtaining a correctly representative 
sample. It was recognized that the distribution of impurities might 
well be affected by factors other than the resistance to crushing and 
that even the crushing factor for the same kind of impurities would 
probably vary with the condition of the coal, for instance with the 
moisture content. Accordingly it was not expected that a law could 
be found connecting the size of particles with the ash content. The 
results of tests tabulated in Table 1 and shown in part graphically 
by figure 1 would seem to be in accord with this reasoning. 

Six samples of coal selected at random were crushed to 20-mesh 
and screened as shown in Table 1, following. Each screened part 
was weighed and analyzed (duplicate determinations) for moisture 
and ash, the ‘‘weighted mean” of all determinations being taken 
as the true value for the original gross sample. The ‘weighted 
mean” was obtained by multiplying the determined percentage of 
moisture or ash by the weight in grams of the sample represented, 
adding the results, and dividing the sum by the total weight of the 
gross sample, The sixth and seventh columns of the table show the 
variations from the true value for cach sized part of the sample. 
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PREPARATION OF SAMPLE OF COAL FOR ANALYSIS. 9 


TABLE 1.—Results of analysis of samples of coal from different stages of sizing. 
BITUMINOUS COAL FROM DUNSMUIR, VANCOUVER ISLAND, B, C. (SAMPLE 57831), 


] 
| Ash in Moisture 
at ; Ash in ; Ash dif- 
Mesh. Weight, , Moisture.! coalas ~ differ- 
t s received. dry coal.| ence. ference. 


Grams. | Per cent. | Per cent. | Per cent. | Per cent. | Per cent. 
272 2.22 11. 98 12, 25 +0, 13 +0,27 


152 2.10 10, 92 nN. 15 + .01 — .&8 
142 2.07 10, 87 11.10 - — .88 
38 1,82 10. 82 11, 02 — .27 — .86 
268 2.00 12.52 12.78 — .09 + .380 
pataerp acces aseecenee'ee [aes saaec% 209" [esseepstcal WHO8 lescesc case |acossacess 


ANTHRACITE BUCKWHEAT (SAMPLE 58077). 


5.30 17.50 18.48 
5.16 17.49 18. 45 
4.81 18.76 19.72 
4.64 i 8 20. 80 
5.01 17 24.40 
5.00 Joccsnccace 20. 64 


ANTHRACITE PEA (SAMPLE 58051), 


4 2.42 

584 2.34 

152 2.73 

102 2.62 

40 2.28 

288 2.17 
Prrrer errs 2.62 


+0.37 +0.55 
— .8 + .11 
— .32 — .39 
— .61 — .4 
— .18 — .47 


PEERLESS NO, 1 POCOHONTAS (SAMPLE 57454). 


OM 40s a totettotetetagtedas 184 0.88 5.61 5.66 | +0.18 | +0.38 

OmG0l.coosascockeecceseacgs tes 108 73 5.00 5.04 | +.03 | —.2%4 

OWS). 2106 oo bokednas caceeus one Ss 63 4.73 4.76 | —.07 | — .52 

ON 100) ccccues sad eucndesecvises 48 17 4.87 4.91 | + .03 | — .37 

| Throigh 1002! S22 saceSszanctes 252 BS 5,32 5.35 | —.12 | + .07 
° Be 


Figure 1 shows these results graphically. The ordinates give for 
each sample, in percentages, the positive and negative variations 
from the weighted mean, which is represented by the horizontal 
base line, and the abscissas show the mesh or fineness of the corre- 
sponding part of sample. 

57956°—18——2 
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10 METHOD OF LEAST SQUARES. 


As regards the anthracite coals, the fines invariably contained the 
most ash and the coarser sizes the least, the ash content of the inter- 
mediate sizes being near the mean value. For the bituminous coals 


Moisture - - - -, ash ——_, 


Figure J.—Graphs of results given in Table 1. 


INGO Udd ‘NOLLVIUVA 


there is no such regularity. Doubtless more data would bring out 
interesting points of regularity for various kinds of coal, but for the 
purpose of this paper the large amount of experimental work necessary 
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PREPARATION OF SAMPLE OF COAL FOR ANALYSIS, 11 


hardly seems warranted, particularly as such data would apply only 
to crushing with rolls. In general, the order of variation that may 
be expected between different sizes is shown, and it is noteworthy 
that the differences are so large, particularly as concerns the coarse 
and fine material of the anthracite coals. In every sample the inter- 
mediate sizes were of lightest weight. 

In this connection it was thought of interest to determine the 
eveness of distribution of the impurities in coal of a given fineness. 
Accordingly two coals were sized, 3-gram samples were taken at 
random from each sized part, and the ash content of each was 
determined. The results, which are given in Table 2 following, 
indicate that different parts of the same mesh sample vary con- 
siderably in ash content. 


TaBLe 2.—Results of tests to determine uniformity of ash content in coal samples of a 
given fineness. 


BITUMINOUS COAL FROM MARION COUNTY, ILL, (SAMPLE 58211). 


| Ash content of grab sample No.— 


Mesh. — a 
| 1 2 3 

Per cent. Percent. Percent, 

9. 57 10, 78 93 

9.25 9.7 8.75 

9.60 9. 9S 10. 65 

9.65 10. 25 10.10 

| 
ANTHRACITE BUCKWHEAT (SAMPLE 68136). 

O40! ts 55.25 26 ert Sesg ashore 21.62 22.38 
On 60 23.45 22.00 
On 80. 23. 88 24.05 
On 100, 24.32 24.65 


For the particular bituminous coal selected the sizing seems to 
have little effect on the ash distribution, whereas with the anthracite 
coal its effect is evident. For the same mesh sample there is, how- 
ever, considerable variation with both coals. 

The importance of maintaining the proper ratio of weight of 
coarsest particles to total weight of sample has been demonstrated 
by Bailey,* who finds by experiment that the relation between 
“size-weight per cent,’’ the probable error, and percentage of 
extraneous ash can be closely approximated by the parabolic equation: 

r=.522,Vap—.353p, 
where r represents the probable error; that is, for a number of repeated trials 
there will be as many errors smaller than r as there will be greater; a represents the 
percentage of extraneous ash, or that coming from impurities in the coal, such as 
slate and pyrite; p represents the size-weight percentage, given by the ratio of the 
weight of the largest piece of slate or bone coal contained in the sample to the total 
Weight of the sample. 


a Bailey, E. G., Accuracy in sampling coal: Jour. Ind. and Eng. Chem., vol. 1, March, 1909, p. 171. 
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This formula serves as a means of determining approximately the 
accuracy with which a sample has been prepared and also as a guide 
in estimating the proper amount of sample to take in any given case. 
Applying the formula to sample 58077 (Table 1), which shows the 
greatest variation of the series due to sizing, the error to be expected 
by the usual method of preparation of the sample may be found as 
follows: (1) Weight of largest piece of slate (20-mesh) equals 0.0036 
gram; weight of smallest sample divided, equals 200 grams; size- 
weight per cent, 0.002, equals p; (2) assume a equals 6 per cent. 

If these values be substituted in the formula, r= 0.522 -/ap —0.353p, 
then r equals 0.056 per cent. This means that by the usual method 
of preparation, an error not greater than 0.056 per cent may reason- 
ably be expected, whereas, if the sample had merely been screened 
with rejection of the coarser parts, an error of 3.76 per cent (Table 1) 
might have resulted. The effect of sizing is therefore manifest and 
the importance of keeping the sample well mixed at all stages of 
crushing and reduction is evident. 

The significance attached to the moisture-content results presented 
in Tables 1 and 2 is that for this method of treatment the moisture 
content varies considerably with the variation in the size of particles. 
One might have expected the fines to hold the most moisture,* 
as none of the samples treated had been freshly crushed. How- 
ever, large variations in moisture content would be caused by 
successive screenings, and hence no conclusions can be drawn. 

It now becomes necessary to discuss briefly the theory of the law 
of error in order that the computations which follow be readily 
understood. A full discussion of the subject will be found in Merri- 
man’s book on least squares.° 

For any given series of observations experience has shown that the 
following statements hold; that they are axiomatic: (1) Small errors 
are more frequent than large ones. (2) Positive and negative errors 
are equally frequent. (3) Very large errors do not occur. 

From these axioms it will be evident that the probability of any 
error happening in a given series is a function of the magnitude of 
the error itself; that is, denoting the error by z and its probability 
by y, the equation y= fz expresses the relation between the two. 
The curve, figure 2, following shows the relation graphically. 

Here y, the probability, is measured along the ordinates and z, 
the error, is measured along the abcissas, the curve being symmetrical 
with respect to the y axis. It will be seen that for every positive 
value for z there is an equal negative value; that when z is small y 
is large, and that when z is very large y approaches zero in accord 


a ¥icldner, A. C., Notes on sampling and analysis of coal: Tech, Paper 76, Bureau of Mines, 1914, p. 23, 
bMerriman, Mansfield, A textbook on the method of least squares, 1901, p. 13. 
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with axioms 1, 2, and 3, given above. fz of the above equation 
admits of determination and may be shown to be equal to 


ke7h?x? (1) 


h and k being constants; k is the value of y when z equals zero and h 
is the measure of precision, so-called. For any series of observations 
where h of the one is larger than h of the other (y being the same in 
both cases), it will be seen that z must be the smaller in the former 
case, This series is therefore the more precise of the two. 


FiGURE 2.—The probability curve, v= J TZ. 


If it is desired to know the probability of an error falling between 
limits z=0 and z=z’, say, it may be found by integration of function 
(1) between limits 0 and 2’, thus: 

Ar 
Sashes eahtdrh (2) 


the constant k being equal to cm Values of this integral will be 
Va 
found tabulated in textbooks on the method of least squares? and will 
usually be found stated as 
2 
= e-" dt where t=hr 
vr 


ut 


0 


Now the probable error, 7, of a series of observations is defined as 
such that “it is an even wager that an error taken at random in the 
series will be greater or less than the probable error.” Its proba- 
bility is therefore one-half. 


o Merriman, Mansfield, A textbook on the method of least squares, 1901, p. 220; Bartlett, Dana P., 
The method of least squares, 1915, appendix, p. 7; Smithsonian Physical Tables, 6th ed., 1914, p. 56. 
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By solving the expression . 
(34 


2 
eae e~dhx 
Vi Jo 
and putting x equal to 7, hr is found to be equal to 0.4769. (3) 


The probability of the occurrence of independent errors 2,, 22, 23, OT 
z, is equal to the product of their separate probabilities, or (from 
equation 1) 


—n 
Pl=h (dr) 9 2 eh? 


from which, if P’ is to be of maximum value, it is required that 22? 
Ud 
be small. Taking the first derivative eas of the expression, placing 


the result equal to zero and reducing, it is found that, n—2h?Z2?=0 
or 


0.4769 
hm als5 5533 but from (3) r= 


whence 


The value z, used heretofore, represents a true error, but the cor- 
responding residual v which represents the difference between the 
arithmetical mean of a series of observations and a single observa- 
tion may not be taken as the true error. The residual v may be said 
to differ from z by some small factor w, = ee av? =v? +u? and v 


oe. sy 
may be taken as approximately equal to — *, from which 27 = 
r, the probable error, is then equal to 
S22 
0.6745. / = 4 
4 ES (4) 


which is the formula used in calculating r from the tables that follow. 

Another value, uv, is also used that may be called the high error for 
a given series. The probability of its occurrence is one in one 
thousand, and it is equal to 4.9 r, r being the probable error, found 
from expression 4. 

To compute the probable error, 7, of a series of observations it is 
only necessary to determine the drithmaticg! mean of the observa- 
tions, find the differences v between individual observations and the 
mean, and then apply formula 4. 

As regards the curves shown in figures 3,4, and 5 (pp. 17, 29, and 32): 

For the curves a, the values of r were those computed from the 
series of observations represented, the h values being found from the 
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relation hr=0.4769 (p. 14). The h values found were substituted in 
the equation, 
y = n—the“¥F 


from which for various assumed values of z the corresponding values 
for y were found and plotted. In order to facilitate plotting, the y 
values found were multiplied by the factors given on the margin of 
the sheets. Where such curves are plotted to the same scale, the 
height of the curves is a measure of their relative precision; that is, 
the higher the curve, the higher the degree of precision represented. 
All the curves shown are comparable except the highest curve in 
figure 5, page 32, which represents the highest degree of precision of the 
three curves and could not be plotted to the same scale as the other 
two. The ordinates shown for these curves equally divide the areas 
included between the curves and the axes, and the z values correspond- 
ing are the probable errors in each case. 

Curves } show the relation between the magnitude of an error z, 
and the per cent of total number of errors for the series falling inside 
limit z. In plotting determined values for ) and r are taken from 
the series represented, whence values of ¢ in the integral 


t 
On} if e-tdt 


may be found for various assumed values of z. Values of P, the 
corresponding probability, may be read from tables of values of the 
integral. Values for P found in the tables multiplied by 100 give 
directly the per cent of total number of errors falling inside the 
limit xz corresponding. The theoretical number of errors falling 
between limits given in Table 4 (p. 19) were calculated by the above 
method. Here the close agreement of the theoretical results with 
those obtained experimentally will be noted. 


ERROR OF THE ASH DETERMINATION. 


The method followed by the authors in making ash determinations 
on individual samples was that recommended by the committee on 
coal analysis of the American Society for Testing Materials,* which 
is as follows: 

Place the porcelain capsules containing the dried sample from the 
moisture determination in a cold muffle furnace, or on the hearth at 
a low temperature, and gradually heat to redness at such a rate as 
to avoid mechanical loss from too rapid expulsion of volatile matter. 
Finish the ignition to constant weight (40.001 gram) at a temper- 


a American Society for Testing Materials, Year Bock, 1915, p. 605. 
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ature between 700° and 750° C. Cool in a desiccator, and weigh as 
soon as cold. 

As regards the determinations represented by the data taken from 
the records (fig. 3), the method mentioned was not followed strictly 
in that the ashing was sometimes started in a hot muffle; further, 
the temperature of the muffle frequently exceeded 750° C. Some of 
the later determinations, however, were made in a muffle provided 
with temperature-control apparatus with which the bureau’s labora- 
tory has recently been equipped. Results reported in Tables 3 and 
4 were obtained by the method of ashing recommended by the com- 
mittee. In this connection, and as a matter of interest, attention is 
called to the limit agreement (0.001), constant weight, as given. An 
analyst using this limit agreement may obtain an error as large as 
the total variation permitted between duplicates if a coal is low in 
ash and contains no carbonates. Further, as is presently shown, the 
limit agreement allowed is not less than the probable error of the 
ash determination as experimentally determined. If it be admitted 
that for certain coals it is difficult to obtain closer agreement, and 
that closer agreement is often unnecessary, it would seem advisable, 
for the sake of consistency, to extend the limit of allowable variation 
to correspond. 

In order to obtain expressions for the accuracy of the ash determi- . 
nation for coals containing less or more than 12 per cent ash, two 
100-gram 60-mesh samples that were known to contain ash as 
required were chosen, and 24 ash determinations were made for 
each. The arithmetical mean of the results was found in each case. 
The residuals, v, were then computed and the probable error found 


by the formula 
N32 
r=.67454/ ae 
n—-1 


as explained on page 14 of this paper. The error u, corresponding 
to a probability of one in one thousand, was also computed, the 
results being given in Table 3 following. 
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18 METHOD OF LEAST SQUARES. 


TaBLe 3.—Results of 24 analyses, establishing ash errors for two 100-gram bottle sam pls 
of coal containing 8.8 and 13.6 per cent of ash. 


o Rejected by Chauvenet's criterion; the deviation from the 7 as first found is r x 3.43, or 0.24 per cent. 


aD 
r=0.67154/°" 0.083 per cent. 


u=4.9r=0.21 per cent. 
Using 7? 
r=0.6745 | xy? . 
V aa O87 per cent. 


u=4.9r—0.33 per cent. 


The results in Table 3 mean that for serics of determinations with 
coals of the same character and ash content as those used in the 
experiments, half of the errors for one might be expected to fall 
within the limit of 0.043 per cent, and for the other within the limit 
0.067 per cent; once in one thousand determinations an error of 0.21 
per cent might be expected for the one and of 0.33 per cent for the 
other. The order of accuracy obtainable is thus defined. The dis- 
tribution of errors for the sample of Table 3 which contains 13.66 per 
cent ash is given in Table 4 (p. 19) with the distribution of errors of 
theseries, 1,000 check samples taken from the records. The theoretical 
number of errors shown in Table 4 were calculated (r, the probable 
error being known) from the formula 


referred to on page 13 of this paper. The values of the integral as 
taken from the tables range from 0 to 1 and each value read for the 
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| | 
: | } 
aetyats | Ash. | Residuale.| Ash. | Residual r).|Residual vy. 
| | 
| 
Per cent. Per cent. Per cent. Per cent, Per cent. 
1 8.92 0.05 13. 70 0.03 0. 
2 8.80 07 13.7 04 
3 8.95 08 13.55 -2 ell ' 
4 8.85 02 13.50 17 16 
5 8.93 -06 13, 60 -07 -06 
6 8.80 07 13. 85 18 19 
7 8, 68 19 13. 75 -08 -09 
8 8.85 02 13. 77 .10 ll 
9 8.90 03 13. 58 -09 
10 8.80 .07 13. 80 13 i4 
ll 8.85 02 13. 82 -15 16 
12 8.85 02 13, 98 4,26 |..... ésvesee 
13 8.87 .00 13. 72 05 -06 
14 8. 93 .06 13. 64 03 02 
15 8.98 et 13. 75 -08 09 
16 8.85 02 13.77 10 ell 
17 8.87 -00 13. 81 04 05 
18 8. 90 03 13, 65 +02 O01 
19 8.93 -06 13.7 - 03 -O4 
20 8.90 -08 13. 60 07 06 
21 8.90 - 03 13.70 -03 04 
22 8.90 . 03 13. 57 .10 09 
23 8.95 - 08 13.72 05 06 
24 8.85 02 13.77 +10 al 
Mean... 8.875 [fBeanae eh Has ehvvaesos 13. 674 13. 66 
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z value corresponding may be taken as a fraction of the total number 
of errors of the series. For example, for the series for which the 
mean ash is 13.66 per cent, the value of h was calculated from equa- 
tion 3 (hr=0.4769) whence the value of hz=t of the integral was 
found. Substituting the value of 0.1 per cent for z (the value of hz=t 
is 0.7 per cent) the corresponding probability is found from the tables 
to be 0.6778. Now there are 24 determinations (and hence 24 errors) 
in the series, and 24 times 0.6778 gives 16 as the theoretical number 
of errors for the series within the limit of 0.1 per cent. Another way 
of stating this result is that 67.78 per cent of the total number of 
errors of the series fall within 0.1 per cent. For the curves in 6, fig. 3, 
the values for the ordinates were calculated for corresponding values 
of a by the method just given. Curves in z, fig. 3, were calculated 


from the formula 
y=7-'he-"* 


as explained on page 15 of this paper, each curve being designated 
by the value of r, the probable error of the series corresponding. 


TaBLe 4.—Results establishing distribution of error in ash determinations. 


[r=0.06 per cent, com puted from a table of values for the probability integral, Pa e-?? dt for argu- 
® id 
ment t= =] 


1000 CHECK SAMPLES WITH ASH UP TO 12.5 PER CENT. 


Limits, per cent. 


OD POH O.1 occ iecins ss see aes 
OG 9 Beas gstasasssacaeiacness 
Oto — .3. 

ho een Sees ieee ee oe, PE 


24 DETERMINATIONS ON 100-GRAM BOTTLE SAMPLE CONTAINING 13.66 PER CENT ASH. 


0.00 to —0.08. 20... pac cecatedctaed | 6 | 8 385 
-00 to — -10... 1.2... ee eee eee 14 16 685 
00 to Nb t55. 2255 oe oajtea seee2 \ 20 19 869 


The significance of these curves may be understood by considering 
that for a given value of z, the number of errors of the series falling 
within limits is proportional to the area included between the curve, 
the axes, and the ordinates for +z and —z. As will be seen, the 
ordinate y measures the precision, or, rather, gives the relative pre- 
cision of the series plotted to the same scale as represented by the 
curves. The steep curves will clearly represent the greater number 
of small errors, and will therefore indicate the greater precision. In 
curves a, figure 3, the curve representing the series r=0.06 per cent 
shows this series to be the most precise of the four represented, prac.- 
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tically all of the errors of the series being included by the curve and 
the ordinates corresponding to values +0.3 per cent and —0.3 per 
cent for the error z, and one-half of the errors within values +0.06 
per cent and —0.06 per cent. The other two curves shown (r=0.067 
per cent and 0.10 per cent) will be seen to represent a lower degree 
of precision. It will be noted that there are two series for which 
r=0.10 per cent and since this is so both are represented by the same 
curve. 

As the foregoing explanation applies to all error-distribution curves 
given in this paper, the theory of such curves, presented subsequently, 
will be clear without further discussion. 


TaBLe 5.—Results establishing error of riffing and of ash determination, 12 successive 
samples of 5-pound 20-mesh sample. 


Ash. 


| Sample. Mean. Moisture. Mean. | Maximum.] Minimum. 


| Determina-| Determina-| 
tion A. | tion B. 
|- : a 
| Percent. | Percent. | Percent.| Percent. | Percent.) Per cent. | Per cent. 
1 6. 55 16.35 | 16.45 | 2.70-2.70 2.70 273 | 2.45 
| 2 16.45 16.55 6.48 | 2.70-2.75 2.73 
3 16. 55 216.95 | 16.75 | 2.75-2.60 2.68 ea 
| 4 16. 80 16.55 | 16.68 | 2.70-2.65 2.68 
| 5 16. 85 16.60 | 16.73 | 2.80-2.70 2.75 
| 6 16. 50 16.55 | 16.53 | 2.65-2.70 2.68 
7 16.50 16.40 | 16.45 | 2.60-2.55 2.58 
8 16.75 16.70 | 16.73 | 2.70-2.55 2.63 
9 16.45 16.40 | 16.43 | 2.55-2.65 2.60 
10 16. 55 16. 65 6.60 | 2. 65-2. 55 2.60 
11 16. 50 16.45 | 16.48 | 2.55-2.50 2.53 
12 16. 80 16.95 | 16.88 | 2.50-2.40 2.45 
}' PMlean 500552 ctl aseee Beeade 16/00. cooetvuctesd eek event sonlsoese 


a Checked by third determination giving 16.75 per ent. 


r=0. 6745/2 =0. 103 per cent. 
(from residuals, Mean = 16.60—mean of determinations A and B for 
samples 1 to 12). 


u=4. 9r=0. 50 per cent. 


Table 5 gives the results of 12 ash determinations of successive 
samples taken from the same 5-pound gross sample, the probable 
error found being 0.103 per cent. This value includes, of course, 
the error of the ash determination and that from sampling. The 
error applies only to the particular coal examined, and in order to 
answer the question as to what the figure would be for the average 
coal containing up to that amount of ash, the records of 1,000 
consecutive check determinations were noted; each of the differences 
between the results for the first sample and those for the duplicate 
sample was halved; and the residuals thus obtained were squared 
and added. The value for r, found by the usual method, was 0.10 
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per cent, or practically the same as that found for the results of 
12 determinations given above. 

It will be appreciated that the 1,000 check determinations do not 
represent 1,000 observations of the same quantity, but rather a 
second observation of each of 1,000 quantities differing slightly in 
value, the mean of each pair of observations being taken as the true 
value. As the errors of observations for all pairs of the series are 
of the same order of magnitude, the law of residuals 


a 
r= 0.67454/ oat 
n—-1 


willapply. In the application of the formula here it should, of course, 
be considered that for every pair of observations there are two 
residuals equal in value and opposite in sign. For example, if one 
observation gives a value of 14.20 per cent and the duplicate 14.40 
per cent, the mean would be 14.30 per cent and the residuals would 
be +0.10 per cent and —0.10 percent. Now if we call n of the 
formula the number of check determinations, and v half the difference 
between duplicates, we must write for the probable error 


je. [2 Te 
r= 0.67454) 92. 


However, if n is large practically the same result is obtained from 


the formula 
[= 
r= sl a, a ad ’ 


which is accordingly used here. 

It might be expected that the result for the 1,000 check samples 
would be lower than that for the individual sample, as the check 
samples included a large number of low-ash coals. However, this 
would be offset to some extent, as variation in moisture content was 
not considered in connection with the 1,000 check samples. The 
curve for this series of 1,000 check samples (r= 0.10 per cent) is shown 
in figure 3, page 17. It will be noted that the series for 12 deter- 
minations (r=.10 per cent) is represented by the same curve. 

The ash in approximately 85 per cent of the coals analyzed in the 
laboratory does not exceed 12 per cent, and it was thought desirable 
to determine the probable error that would represent coals containing 
not to exceed 12 per cent ash. The figure was found to be 0.067 
per cent calculated from the records of 1,000 check determinations, 
the method of computation being the same as that for the 1,000 
check determinations referred to above. However, only determina- 
tions for the same bottle sample were considered in this calculation. 
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Data for 106 sets of duplicate samples analyzed in the bureau’s 
laboratory at Pittsburgh, Pa., are given in Table 6, following.¢ 

The differences between duplicates given in this table were halved 
and the results squared. The sum of these squared residuals was 
found to be 151.2650, which may be taken as 2? of the formula 
for the probable error r, 


Substituting the value found and 105 for n—-1 of the formula and 
solving, r is found to equal 0.083 per cent. The value of the error 
corresponding to a probability of one in one thousand (4.9 r) is 
found to be 0.406 per cent. 


TaBLe 6.—Results establishing the ash error in the analyses of 106 sets of duplicate 
samples. 


Ash in moisture-free coal. 
Labora- Kind of coal, Source. 
tory No. 14 Reserve : 
Original. f Difference. 
s I 

Per cent. Per cent. Per cent. 
10900 | 6. 99 7.30 0. 31 
11066 |... 9. SA 9,90 36 
1lw9l4 9. 24 9.33 09 
11361 | 7.91 8.01 10 
11294 11.15 10. 94 21 
11048 | Subbituminous......... Converse County, Wyo.......- 6. 42 6.35 07 
11059 Joo... WMO crcwtcccdcadseetaa | Johnson County, Wyo........ 18. 02 18. 32 30 
10827 |.....do.. } 1 6.75 7.18 43 
LWT 1 | 7.60 7.70 10 
10825 |.....d9.. 5.88 6, 22 34 
12006 4.70 4.86 16 
Uinta County, U 2 9.17 14 
Wasatch County, Utah. 3.5) 3.7 13 
Oye teow eh . 6 7.70 02 

El Pas 20.7 21.33 : 

Ny nae do 8.75 8.71 

Philippine Islands 8.17 8.69 52 
Hot Springs County, Wyo 6.17 6.08 14 
a (7 . 5.94 oS eee 
5. 80 5.7, 05 
5.43 5.54 ul 
5.72 5. 46 26 
5. 92 5. 66 26 
5.70 5. 89 19 
.| Jefferson County, Ala. 10. 91 11.38 A? 
Las Animas County, C | 11.75 11. 92 17 
Montgomery County, Il 11.47 11.45 +02 
Saline County, Hl... 2... 7.98 7.77 +21 
Vermillion County, Il... 12.02 12. 29 +27 
1G cs AOS ip betath tsb ouck 11.3 11.95 +12 
waste aceaehace tee super enens, 12.07 12.04 , 8 
Pike County, Ky...2.....0..0. 2.75 3.04 29 
Henry County, Mo.. 16. 80 16. 78 02 
pS: resseeacadeseecesbereene 14. 26 14. 23 +03 
setUOtet isch. ssapeeecacesataie 11. 66 11. 96 +30 
Randolph County, Mo.......- 9. 28 9.19 09 
Pittsburgh Couniy, Okla...... 4.57 4.24 33 


a wielanet, A.C., Notes on the sampling and analysis of coal: Tech. Paper 76, Bureau of Mines, 1914, 
pp. 53-35. 
‘ + The sample received at the laboratory, after having been crushed to pass a 10-mesh sereen, is aie 
nto two parts, one part being pulverized in the ball mill fur the regular sample, and the other t 
placed back in the can and held in reserve. -: nse i 
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TaBLE 6.—Results establishing the ash error in the analyses of 106 sets of duplicate 
sam ples—Continued. 


Ash in moisture-free coal. 


Labora- Kind of coal. ; Souree. 


tory No. Original. ry Difference. 


Per cent. Per cent. Per cent. 


12. 49 12.90 0.41 
15. 57 15.41 -16 
11. 94 12. 02 08 
11. 96 12. 54 - 58 
12. 50 12. 80 30 
11.95 12. 08 -10 
10, 24 10. 40 -16 
10. 29 10. 43 14 
9. 80 9. 80 00 
6.13 5.99 -14 
5.92 5.97 +05 
9. 21 9. 30 -09 
21.31 21. 23 - 08 
11. 24 11.12 12 
11.13 11, 25 12 
10. 99 10. 86 13 
10. 93 10. 98 - 05 
11. 04 1 Pe Borniccecce 
11.00 BE OO asc cadiecese 
11.08 11.01 07 
11. 23 11.16 07 
7. 52 7.47 05 
6.89 7.25 . 36 
idhea County, Tenn 12, 03 11.98 05 
Carbon County, Utah......... 7.37 7.04 33 
Russell County, Va..........- 5. 93 5.81 -12 
aise COs. 2 ecce cou wbaweriestinakseee 5. 44 4.86 58 
King County, Wash........... 11.38 11. 54 +16 
| Pierce County, Wash.......... 15. 32 15. 00 32 
Weston County, Wyo.......-- 16, 84 17.10 26 
263 Denacestut Siaca sec cscmasce 18, 40 18. 60 20 
West Virginia 7.94 7.72 22 
rele do..... 6.78 6. 60 18 
Marion County, W. 6. 08 6. 25 17 
| Marshall County, W. Va | 7.26 7.16 10 
Japan | 5.03 4.75 28 
.-do 5. 683 5. 69 06 
--do. 13. 37 13, 01 36 
passceopocsar ies do 13. 20 13. 30 -10 
OWING. i925 e025 cececsciezsncs ¢ 15. 92 16. 39 47 
Westmoreland County, Va.... 9. 58 9. 47 “u 
Cambria County, Pa.........- 6.15 6, 45 30 
New aiver, Fayette County, 2.83 2. 56 27 
Ww. 
New River, Raleigh County, 3.65 3.58 07 
W. Va. 
eine (| ee reece 3.7 3. 64 06 
Pocahontas “bed, McDowell 11.32 11.38 - 06 
County, W. Va. 
Penns yi Vans « ss ssc es esnsssets | 12, 22 12. 03 -19 
di 25. 45 25, 08 42 
7. 23 7.05 18 
6.03 5.89 4 
11.94 “11.66 28 
13952 | Anthracite, No. 1 buck- 17.55 17.45 10 
17.11 37 
17.27 24 
17, 21 .19 
17.35 25 
17.53 7) 
18. 11 -73 
17.33 - 03 
13905 | Anthracite, No. 3 btick- 18. 37 28 
wheat. 
18. 69 .02 
19.00 1 
18. 84 12 
18, 43 .32 
18. 56 -14 
17.49 -00 
17.48 . 03 
17. 60 . 02 
17.96 «US 
\ 
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This series, as will be seen, showed slightly higher precision than 
was shown by the series of 1,000 check determinations including 
coals containing up to 17 per cent ash; however, in the latter the 
moisture variation was not considered and consequently the result 
might be expected to be somewhat high. 

Now in this series, for which r equals 0.083 per cent, there were 
five pairs of duplicates where the variation was as much as 0.5 per cent. 
It will be interesting to calculate for the series the number required 
by theory. A variation of 0.5 per cent between duplicates means 
that for each determination the variation from the mean is 0.25 per 
cent; therefore we may consider that we have 10 determinations where 
the error is 0.25 per cent. Knowing the value of r (0.083 per cent) 
for the series we may read from the tables of the probability integral 


t 
“a f e~ dt, any value of x of the series. In this instance, where z 
To 


equals 0.25 per cent, the probability is 0.9577, which means that 95.77 
per cent of the errors of the series are within limit 0.25 percent. There 
are 106 duplicate determinations for the series, hence 212 single de- 
terminations and 212 errors. As 95.77 per cent of 212 equals 203, 
203 errors less than 0.25 per cent are required for the series by 
theory; that is, 101.5 pairs of duplicates. Then 106—101.5=4.5 
which should vary more than 0.5 per cent, required by theory, whereas 
5 were actually found. Clearly the agreement is close. A wider 
difference between theory and practice might have been expected 
because the average mean used must have differed somewhat from 
the true mean and because the true ash can not be accurately defined. 

It would seem almost axiomatic that determinations for coals of 
high ash content would be subject to greater error than those for 
coals of low ash content, although it is easily seen that the ash dis- 
tribution, and therefore the difficulty with some coals of obtaining 
a sample representative of the impurities present, might alter this 
conclusion. That the proposition is true in general, however, is in- 
dicated by the following series: 

(1) For ash under 5 per cent, a series of 42 check determinations 
gives r equals 0.027 per cent. 

(2) For ash from 5 per cent to 10 per cent, a series of 103 check 
determinations gives r equals 0.033 per cent. 

(3) For ash from 10 per cent to 15 per cent, a series of 66 check 
determinations gives r equals 0.124 per cent. 

(4) For a series of 67 check samples in which the ash was 15 to 25 
per cent the probable error was 0.20 per cent. 

Series 4 includes the error of the ash determination and that of 
sampling, and may be taken to indicate about what may be expected 
for check samples of coals containing the maximum amount of ash. 
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The method followed in the Bureau of Mines laboratory in making 
& moisture determination is as follows: 

One gram (+0.2 milligram) of the sample is weighed into a No. 2 
porcelain capsule which is provided with a numbered thin aluminum 
lid. The capsule, uncovered, is placed in a glycerin-jacketed oven 
and dried at a temperature of 105° to 110°C. for 1 hour. During the 
drying, warm air, previously dried by bubbling through concentrated 
sulphuric acid, is passed through the drying compartment of the 
oven at the rate of about 4 volumes a minute. The capsule contain- 
ing the dried coal is then taken from the oven, covered, cooled in a 
desiccator, and weighed. The first weighing of the sample, as well 
as the transfer and weighing of the dried sample, is effected as quickly 
as possible in order to prevent the sample from taking up moisture 
through exposure. In this connection the authors wish to state that 
they have found it comparatively easy to weigh to 0.2 milligram 
without introducing an error of greater magnitude through changes 
in moisture content of the sample. 

The weight of the sample as well as subsequent weights may be 
quickly obtained by ‘‘snap weighing” or weighing to a standstill— 
that is, adjusting the weights till the balance pointer remains at 
zero without swinging. 

Reasonable accuracy in “snap weighing’”’ requires considerable 
practice; also a quick accurate balance is absolutely essential; its 
sensibility should be high. Most analytical balances are not satis- 
factory for “snap weighing”; they are too slow and do not permit 
the required delicate adjustment. A light balance that will weigh 
50 to 100 grams and is sensitive to one-twentieth milligram is to be 
preferred and it must be so adjusted that it will quickly indicate an 
excess of 0.1 milligram from standstill. It is well, of course, to prac- 
tice making ‘‘snap weighings’”’ with some substance not subject to 
moisture changes, as the weight first obtained may then be checked. 

Table 7 following shows the results of two series of determinations 
for moisture—one series for a coal having a moisture content of 1.03 
per cent and one for a coal containing 8.50 per cent moisture. 
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TaBLE 7.—Results establishing the moisture error for two 100-gram bottle samples 
containing 1.03 and 8.50 per cent moisture. 


| 
Series 1. | Series 2. 
Test No. Se 
Moisture. | Residual v.| Moisture. | Residual v. 
| | 
| 
| 0.98 —0.05 8.55 0.05 
1.02 — .01 8.59 09 
92 -.1l 8.45 — .05 
98 — .05 8. 50 00 
95 - OS 8.45 — .06 
1.00 — .03 8.55 -05 
1.93 00 8.60 .10 
98 — .05 8.61 ll 
1.06 203 8.65 215 
115 212 8. 60. -10 
119 07 8.43 — .07 
1,06 203 8.40 — .10 
1.10 07 8.48 — .02 
1.05 «02 8.45 — .05 
1.05 «02 8.33 -. 
1.05 02 8.35 — .15 
1.05 02 8.0 10 
1.95 02 8.58 . 08 
1.07 - 04 8.60 .10 
98 — .05 8.45 — .05 
1:02 — .ol 8.38 -.R 
1.07 04 8.45 — .05 
1,05 02 8.40 — .10 
TOGA secetaedssex | B S02}. oi seksses es 
\ | | 
Yet 
r=0.67454 /——-=0.034 per cent. 
n—-1 
u=4,9r==0.17 per cent. u=4.9r=0.31 per cent. 


The distribution of errors for the series of tests covering the sam- 
ple with a moisture content of 8.502 per cent was computed in the 
same manner as the distribution of errors shown in Table 4 (p. 19), 
being based on the following values: 


Results establishing distribution of error in analyses of 100-gram bottle sample of cocl 
containing 8.502 per cent moisture. 


Theoretica! 
are : Errors | Theoretical 
Limits, per cent. found. muumber, wrrood in 
OM COO OF as 5 ae gaa taadanac’ascens 8 | & 385 
10 ; 7 685 
05 to— 15 | a) 19 | Soo 


All the determinations were made from the same 100-gram bottle 
sample as was used in the ash determination represented in Table 3 
(p. 18), and there seems to be nothing deserving of particular motice 
except that there is no preponderance of negative residuals with the 
last determinations of either series. It might be expected that the 
last samples weighed would run low in moisture, as it was necessary 
to expose the whole sample repeatedly while the first samples were 
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being weighed; and it may be stated further that none of the samples 
had been previously air dried. It would seem, then, that there is in 
reality not much moisture change in the gross sample during the 
weighing of the moisture sample. 

It is the practice of the Bureau of Mines laboratory to make a 
separate moisture determination from a 20-mesh grab sample taken 
immediately after the 5-pound sample has been crushed to 20 mesh, 
a 5-gram sample being weighed out for the determination. No figures 
were obtained for the possible precision with this 5-gram 20-mesh 
sample, but it was noted in a general way, while the data on 400 check 
determinations for the 60-mesh sample were being taken from the 
records, that the agreement for the determinations with the 20-mesh 
sample was better than that for the 60-mesh sample. As the 20-mesh 
sample almost invariably contains the most moisture, some of which 
may be extraneous,* it would seem that, in general, errors of the 
moisture determination are made up mostly of errors in weighing 
and are not largely due to gain or loss through exposure. The 
weighing error connected with the 5-gram 20-mesh sample is, of 
course, divided by 5, whereas an exposure error might be expected 
to vary as the area exposed. 

If the size of the particles can be said to have an influence on the 
tendency of a coal to retain or lose moisture during exposure, the 
20-mesh sample might be expected to change more rapidly than 
the 60-mesh sample. As in reality the 20-mesh samples considered 
here consist mostly of finer coal (Table 1, p. 9) the effect of the coarser 
particles would be to keep the sample from packing readily and thus 
to leave more of the fine particles, which are conceived to retain the 
most extraneous moisture, exposed to air circulation. 

The probable error for the 60-mesh coal in 400 check determina- 
tions represented in the records was 0.034 per cent, 95 per cent of the 
determinations showing a moisture content of less than 4 per cent. 
This degree of precision, as will be seen, is about the same as that 
obtainable in the ash determination with coals containing 5 to 10 
per cent ash. Seemingly the crror for the moisture determination 
increases more rapidly with increasing moisture content than does 
the ash error for a corresponding increase in ash content. 


ERROR IN VOLATILE-MATTER DETERMINATION. 


The method used by the authors for the determination of volatile 
matter was as follows: 

One gram of the 60-mesh coal was weighed into a 15-c. ¢. platinum 
crucible provided with a closely fitting cover and heated for seven 
minutes in a small electrically heated muffle furnace kept at a tem- 


®Fieldner, A. C., Notes on the sampling and analysis of coal: Tech. Paper 76, Bureau of Mines, 1914, p.23, 
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perature of 950° C. (415° C.) by means of an automatic temperaure 
regulator. The crucible containing the residue was then cooled on a 
heavy cast-iron plate and weighed. The loss in weight minus the 
moisture found for the sample was taken as the volatile matter. The 
crucibles were always placed in the same relative position in the 
furnace, being placed as near as possible to the regulating thermo- 
couple. However, the temperature in a considerable part of the 
furnace does not vary appreciably from that at the point of the 
regulating couple, as the authors have demonstrated by exploring 
the furnace with a quick-acting auxiliary pyrometer. 

The results of a series of 24 determinations on a 100-gram 60-mesh 
sample of semibutuminous coal are given in Table 8 following: 


TaBLeE 8.—Results establishing volatile-matter error for 100-gram bottle sample of 
semibituminous coal having a volatile-matter content of 19.32 per cent. 


| Volatile Volatile | 
Test No. matter, | Residual v. || Test No. matter, | liesidual y. 
| per cent. | per cent. 
>. | i 
19.75 0. 43 19. 80 0.52 
19. 50 18 19. 30 —.02 
19.00 —.32 19.15 -.17 
19. 60 19.55 -B 
19.70 38 19,10 —.22 
18. 95 —.37 19.75 43 
19.30 —.02 19. 25 —.07 
19. 25 —.07 19. 25 -.07 
19.65 33 19.10 —.02 
19.45 13 19.30 —.02 
19. 50 18 19.20 —.12 
19.30 —.02 19.32 
19.15 -.17 
ry 
r=0.6745 n—17—0-18 per cent 
u=4.9r=0,86 per cent, 
Distribution of errors, 
. ; Theoretical} 
Errors | Theoretical! ; 
Limit, per cent. found. | number. nabs an 
’ 
0.10 to —0. 10... 
»20to— .W.. 
.30to— .30.. 
40 to — .40 


Curves based on the results presented in Table 8 are shown in 
figure 4. A series of 500 check determinations contained in the 
laboratory records gave a probable error of 0.104 per cent. That 
the error for the large series is less than that for the 24 determina- 
tions on the same sample is perhaps due to the fact that the bitu- 
minous coal used in the smaller series gave off its volatile more 
rapidly than does the average coal, causing greater variations. In 
this connection it may be stated that for the series of 500 check 
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determinations no determinations originally rejected were considered. 
The determinations rejected represented tests in which explosions 
occurred or in which other errors were evident. The effect of an ex- 
plosion in the crucible is to lift the lid of the crucible, thereby admit- 
ting air which oxidizes somewhat the coke residue, making the result 
for the determination too high. 

In the experience of the authors the use of the electric furnace is 
attended with more trouble from ‘‘explosions” than is the use of the 
Bunsen burner. With the furnace the temperature within the 
crucible comes to a maximum more quickly, and any large amount of 
volatile matter present is evolved with explosive force, whereas with 
the burner the evolution of volatile matter takes place more slowly. 
All things considered, however, the furnace gives the best results as 
it is capable of closer temperature regulation. For constancy of 
results it has been shown @ that control of the temperature of the 
distillation within narrow limits is essential, and although a method 
might be evolved that would obviate the explosion trouble for a par- 
ticular kind of coal, such a method would probably be subject to more 
serious error in other respects. 


ERROR IN HEATING-VALUE DETERMINATION. 


In the Bureau of Mines laboratory the heating-value determina- 
tion is made with an acid-proof oxygen bomb and a water-jacketed 
calorimeter designed especially for the bureau’s work. The pre- 
cision attainable is somewhat higher than that obtained with the usual 
form of fuel calorimeter; further, the time required for manipulation 
and observation is considerably shorter than with many other types. 

Owing to the seeming simplicity of the heating-value determina- 
tion, the results obtained in the determination may be assumed to 
be of unquestioned accuracy when such an assumption is not war- 
ranted. The determination requires fully as much manipulative 
skill and close attention to details as does the carbon-hydrogen de- 
termination, and the effect of the personal equation is therefore apt 
to be large in comparison with the effect of errors attributable to the 
instrument used. 

The results given in Table 9 following are regarded as indicating 
the best work possible with the bureau’s instrument in the hands 
of a particularly skillful manipulator. 


a Fieldner, A. C., and Hall, A. E., Influence of temperature on the determination of volatile matter in 
coal: Proc. Eighth Internat. Cong. App. Chem.,, vol. 10, 1912, p. 139, 
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TaBLeE 9.—Results establishing the heating-value error in 24 determinations on boitle 
sample of 100 grams of coal, 


Date of determinat ion. 


Heating Residual y2 
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value. 
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r0.07854/; — =3.2 B. t.u. 
u=4.9r=15.7 B. t. u. 


The curve representing the results presented in the table is 
given in figure 5. There are also plotted in figure 4 the curve for 
1,000 check determinations where r=16 B. t. u. and the curve for 
500 check determinations where r=24.7 B. t. u. The curve r=16 
B. t. u. includes variations due to the personal factor, and to differences 
in calibration of the different instruments, and indicates the precision 
to be-expected under the conditions of regular routine work. The 
curve for r=24.7 B. t. u. includes all errors considered in the series 
7r=16 and in addition the widest variations occurring in the series, 
no determinations being rejected; that is, a number of observations 
known to have large constant errors were included. This series 
indicates merely what would happen if no rule for the rejection of 
doubtful observations were available and no attention were paid to 
such sources of error as slight amounts of unburned coal, thermostat 
setting slightly in error, excessive contact resistance in the ther- 
mometer switch, or a mistake of 0.001 ohm in the bridge reading. 

As the determined heating value is a function of the moisture and 
ash content of the sample, it is of interest to calculate the heating 
value of the coal substance and to determine the effect of the combined 
errors for ash, moisture, and heating value on the calculated result. 
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If variations from the true ash and moisture values be neglected 

the heating value of the coal substance is given by the formula: 

a 

100—(z+y) 

Where B=calculated heating value, on ash-free and moisture-free basia 
A=heating value as determined. 


=value found for ash. 
y=value found for moisture. 


(1) B= 


The following assumption may be made: 


r,=0.15 per cent for the probable error in A. 
r,=0.10 per cent for the probable error in x (including sampling error). 
r,;=0.04 per cent for the probable error in y (including sampling error). 
R=the combined errors. 

The combined error f, corresponding to the error in B, is given by 


the following formula, @ 


AES 5 Ne Ae Oe 
(2) r=n[ 54 | +r] 5 | +r3| 5 


applied to the formula given above. 
Differentiating: 


1 2 A 2 2 
« 2» 2 = a 2 = a +) es “= 
(3) ae [ W222 5 | tT, [omega +s [; ma 


Now, assuming A= 14,000, x =8.00 per cent, and = 2.00 per cent, 
and substituting these values and the values found for r,, r,, and r, 
in equation 3, we have &=-/0.00000272 + 0.00477 +0.0298 =0.18 
(nearly), which is the combined error for a coal of the character 
assumed. Of course the result would vary for different coals, but 
the foregoing calculation shows that the error introduced by the 
combination is of the same order as that of the individual errors 
and is not at all serious. Variations found, then, for the calculated 
heating value of the coal substance must be due largely to variations 
in percentages of volatile inorganic constituents (moisture, CO, and 
SO,) of the ash,’ variations in composition of the coal substance it- 
self, and, perhaps also to some degree, a variable affinity for moisture. 

It has been shown‘ that the heating value of coal substance as 
calculated from analyses of coal samples from the same seam is 
subject to variations usually not greater than 1 per cent.? 


CONCLUSION. 


The results obtained through the application of the theory of proba- 
bility would seem to indicate rather clearly what degree of precision 
may be expected with the coal analysis when the quantities determined 


a Bartlett, D. P., The method of least squares, 3rd ed., 1915, p. 57; Merriman, Mansfield, A textbook on 


the method of least squares, 1901, p. 79. 
> Parr, 8. W., Chemical composition of Illinois coals: Illinois State Geol. Survey Bull. 16, 1910, p. 22. 
¢ Fieldner, A. ‘c. ., Notes on the sampling and analysis of coal: Tech, Paper 76, Bureau of Mines, 1914, p. 47, 
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are defined by the actual methods used for the analysis. When the 
error corresponding to a probability of one-half is taken, the errors in- 
volved in sampling and analyzing coal may be taken to be as follows: 
In taking the 5-pound gross sample, about 0.04 per cent of sample, 
because of failure to include or exclude impurities represented, is 
ash in the analysis; in the ash determination, 0.08 per cent; in the 
moisture determination, 0.04 per cent; in the volatile-matter deter- 
mination, 0.15 per cent, and in the heating-value determination, 
0.15 per cent. Combining these possible errors makes a possible 
error in calculating the heating value of coal substance of about 
0.18 per cent of the heating value. Now the limit determining 
rejection is often taken as the error whose probability is 1 in 1000, 
which is 4.9 times (or roughly 5 times) the probable error of a single 
determination; the outside limit for error in the above determinations 
may be easily calculated. The limits of error would then be as follows: 
For sampling, 0.20 per cent; for the ash determination, 0.40 per cent; 
for the moisture determination, 0.20 per cent; and for the heating- 
value determination, 0.75 per cent. It would seem, then, that. the 
limits allowed by the committee on coal analysis of the American 
Society for Testing Materials (see p. 6) are not limits outside which, 
with ordinary care, determinations would never fall; it would seem 
rather that the committee means to designate limits within which a 
large percentage of the errors will fall. 
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